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>>. 1 Introduction 

Localized surface plasmons (LSP) of metal nanoparticles [T], [2], [3], |1] enjoy a potential 
as elements of nanooptical networks. These resonant modes arise from the excitation 
\ of coherent oscillations of conduction-band electrons, localized on the surface of these 
^ ! particles. The strong light coupling leads to its absorption and spatial confinement to a 
CN ■ nanometric scale, which results in large local enhancements of electromagnetic field in- 
^ ■ tensities. Ensembles of particle chains have been extensively studied [S] . That provides 
nanoscale control of the transmission, manipulation and switching of optical signals [6] . In 
addition, nanoparticles either dielectric [7] or metallic |8] have been studied in microdisks 
in connection with the mutual perturbation of their resonances, such as counterpropa- 
gating whispering gallery mode (WGM) splitting, modification and their monitoring in 
the cavity, or control of the radiating properties of the nanoparticle by the WGMs of the 
cavity [9]. 

A somewhat related phenomenon known as enhanced optical transmission through sub- 
wavelength apertures, either alone or forming a grating, [10], [TT], [12], [13], [H], has 
received much attention in connection with its potential application for light concentra- 
tion, detection and wavefront steering. 

These morphology dependent resonances (MDR) also play an important role in photonic 
crystals (PC). For example, their influence on the crystal bandgap size and position has 
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been studied [15], [16]. In this way, the Mie resonances of the particles forming the PC 
constitute the hght propagation and enhancement vehicle in the upper bands of those so 
called molecular photonic crystals [T7]. On the other hand, as regards metallic PCs, there 
are certain advantages in fabricating them, like reduced size and weight, easier production 
methods and lower costs, as well as the fact that low loss metal PCs have been studied 
[T8] . Furthermore, many applications are being developed from these structures, like 
waveguiding [TH], [20], [21], [22], waveguide mode-plasmon coupling [23], light transmission 
control [21], thermovoltaics and blackbody emission [25], [2S], and lensing [27]. 

In this paper we wish to answer two questions: first, what is the effect of the presence of 
plasmonic particles in the transmission zone of a subwavelength aperture which is super- 
transmitting? Do these particles enhance or inhibit this supertransmission?. Second, how 
does the energy flow of transmitted light by this slit, (or slits), propagate through these 
sets of particles?. We present a study by means of numerical simulations that show new 
effects in configurations of plasmonic nanoparticles near the exit of supertransmitting sub- 
wavelength slits practiced in a thick slab. Our study is carried out in 2D but the essential 
features observed as regards enhanced transmission and coupling of resonances and light 
transport are likewise obtained in 3D [21] [29]. Also, this 2D geometry constitutes a good 
model with equivalent effective constitutive parameters for microdisks [30] , [31] • Further- 
more, such a geometry is adequate to deal with structures of long parallel nanocylinders 
in 2D PCs or metamaterials [32]. We then address localized surface plasmons of metal- 
lic nanocylinders [33]. We see the behavior of light concentration and transmitted field 
enhancement efficiency in this kind of particle sets (which appears maximized when the 
configuration of the plasmonic set has properties similar to those of a nanoantenna) . 

Further calculations deal with configurations like linear and bifurcated nanocylinder 
chains, also addressing the natural step passing from particle chains, placed one close 
to another, towards a photonic crystal (PC) geometry situated in front of a metallic array 
of slits. In this way, we design a method of collimation and coupling of light from free 
space into the particles by placing them close to the slits. These numerical simulations 
allow to study the effects that arise in the near field as regards extraordinary transmis- 
sion in the slits and the excitation of LSPs of these metallic particles; and since the 
calculations are exact, they constitute a rehable design of future experiments that can 
be performed in either 2D or 3D. In this respect, our calculations indicate that whereas 
the nanoaperture behaves as a collimator that transmits light into the nanoparticles, the 
fundamental role as regards enhancement and field concentration, corresponds to the ex- 
citation of the particle LSP which couple with the MDRs of the slits. This is thoroughly 
studied by first displaying the simple configuration of a single plasmonic nanoparticle in 
front of a nanoaperture, and afterwards extending these observations to more complex 
sets of these particles, including chains and PCs. We take advantage, both in chains and 
in PC geometries, of the excitation of dipolar LSPs. 
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2 Transmission into plasmonic nanoparticles through a nanosUt 
2.1 Numerical procedures. 

From now on, all refractive indexes under the different wavelengths on use are taken 
from |34] and [35]. All particles in this study are considered of Silver (Ag), (refractive 
index n = 0.188 + il.610 at A = 366nm and n = 0.233 + 21.27 at A = 346nm) be- 
cause of their rich LSP spectrum in the near ultraviolet. However, it should be stressed 
here that this is done for the sake of illustrating the effects, and that other noble metal 
material can be chosen. Since the slab is thicker than usual in supertransmission exper- 
iments, in order that the slit walls present high reflection and as small as possible skin 
depth and losses, the metal of the slab is assumed to be Tungsten (W) (refractive index 
n = 3.40 + i2.65 at A = 366nm and n = 3.15 + i2.68 at A = 346nm). It should be re- 
marked in this connection that, ideally, a quasi-perfect conducting slab would exhibit the 
most pronounced effects under study, but if experiments are done with thinner slabs these 
other materials like Al, Au or Ag may be employed. The conservation of light polariza- 
tion is ensured by launching on the aperture linearly polarized light beams of rectangular 
profile, their widths being that of the simulation window (this one always coincides to 
the slab width). The direction of propagation of such beams is normally incident to the 
axis of the infinite cylinders, and their sense is, in all cases, that from the bottom to 
the top of the window. The wavelengths utilized to visually show the slits as well as 
the particle - slit configurations have been chosen so as to maximize their light transmis- 
sion. The distances between all of such particle structures and the plane exit of the slit, 
optimizing their responses, have been expressed in semi-multiples of the corresponding 
LSP resonance - to - supertransmission matched wavelength (with the wavelength reso- 
lution limited by [35]). Maxwell equations are solved by using a finite element method 
(FE) (FEMLAB of COMSOL, http://www.comsol.com) [36], [37]. The solution domain 
is meshed with element growth rate: 1.55, meshing curvature factor: 0.65, approximately; 
the geometrical resolution parameters consist of 25 points per boundary segment to take 
into account curved geometries so as to adapt the finite elements to the geometry and 
optimize the convergence of the solution. The final mesh contains about 10^ elements. To 
solve Helmholtz equation, the UMFPACK direct is employed, the results being displayed 
in stationary regime of propagation. The boundary conditions of the simulation space 
are properly set both to keep the calculations from undesired window reflections and to 
avoid possible geometrical discontinuities. We hence ensure that neither inconsistencies 
due to properties discontinuities of the objets under study nor possible systematic errors 
because of the simulation window interfere with the field calculation. We always select p 
- polarized incident waves to seek extraordinary transmission in this 2D slit (s - polarized 
waves do not produce such a phenomenon in 2D subwavelength slits) [12], [H], [38]. The 
results are thus expressed in terms of either the magnetic vector H(r)[/l/m] which is 
along the cylinder OZ - axis, the electric vector E(r)[l^/m] and the time average energy 
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flow < S(r) > [J/(m^ ■ s)], these last two being both transversal, namely, in the plane of 
the images to show next. The incident wave is a tapered rectangular plain wave of unit 
amplitude in |H| = lA/m{ST), this value corresponding to < S >~ 190W / {m^) . Finally, 
the nomenclature followed, further on, to classify the LSP resonances of the cylinders will 
use the subscripts (i, j), i and j standing for their angular i - th and radial j - th orders, 
respectively. 

2.2 One metallic particle in front of a nanoslit 




Figure 1: (a) Magnetic field modulus |H(r)| (in A/m units), in a W slab aperture (re- 
fractive index n = 3.39 + i2.66, slab width D = 2850nm, slab thickness h = 259.40nm, 
slit width d = 43.23nm). The magnetic vector of the incident radiation (A = 364. 7nm) is 
p - polarized, of unit amplitude and impinges on the slab from below; (b) Time average 
energy flow < S(r) > (in J/{m? ■ s units), maximum arrow length '^2>.'&2KeV / {nm? ■ s), 
minimum arrow length ^ {]eV/{nm^ ■ s)), both norm (colors) and directions (arrows) are 
shown under the same conditions as in Fig. 1(a); (c) Detail of the saddle point evidenced 
by < S(r) > inside the slit (the arrows here appear normalized to their magnitude which 
is shown according to the color bar); (d) Detail of the electric field E(r) (in V/m units), 
both its norm (colors) and directions (arrows) are shown under the same conditions as in 
Fig. 1(a). 
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The phenomenon of field transmission enhancement through a subwavelength aperture 
by excitation of LSP resonances in nearby plasmonic particles, is studied by comparing 
it with the transmission through a slit alone, practiced in a metallic slab. We chose 
A = 364. Tram for the incident wave, where the aperture is supertransmitting. Figure 1(a) 
as well as Fig. 1(b) show the magnetic field norm |H(r)| and the time averaged energy 
fiow < S(r) > distributions in such a slit, respectively (in all figures, from now on, the p 
- polarized wave illumination incides upwards). These distributions, which correspond to 
an incident wavelength at which the slit produces extraordinary transmission, show two 
interesting features: first, the waveguide p - eigenmode inside the aperture and, second, 
the change of direction of the energy fiow, from the lower region of the aperture at which 
the energy is partly refiected, to the upper region of the slit where the energy manifests 
transmission upwards. This sign change is evidenced in Fig. 1(c) by the energy fiow as a 
potential saddle point due to a change of sign of the magnetic field H(r). On the other 
hand. Fig. 1(d) exhibits the electric field E(r) which shows a strong charge concentration, 
and a resulting a dipolar pattern configuration at the corners of both the entrance and 
the exit of the slit (see the two upper and lower vertices of the aperture in this figure). 
Incidentally, we believe that this high intensity concentration at the edges of the aperture 
exit is responsible for the existence of gradient optical forces on dielectric nanoparticles 
placed in its proximity, and the corresponding creation of two potential wells, each of them 
in front of each corner of the aperture exit [cf. Fig. 1(d)], as observed in [32] • Notice in 
those figures that the supertransmitted intensity by the slit is already seen in the spatial 
distribution of |H(r)|, E(r) and | < S(r) > | in the neighborhood of the exit of the slit, 
(we recall that incident |H| = lA/m, |E| = 380V/m and | < S > | = lQQW/{w?). 

As a first example of the effect placing a plasmonic particle near the aperture exit. Fig. 
2(a) deals with an Ag nanocylinder placed in front of the exit of a subwavelength slit 
in a W slab, which presents a large transmission at the chosen wavelength. As shown, 
at A = 366nm and at a distance between the cylinder and the slab d = 549nm, a 
strong intensity enhancement of |H(r)| appears on the cylinder surface, manifested by 
the LSP stationary interference pattern surrounding it, and a standing wave pattern due 
to refiections between the slab/aperture and the particle is observed. It should be noticed 
that this cylinder alone has this LSP^i mode at A = 359. 7nm which is red - shifted as 
seen in Figs. 2(a), 2(c) and 2(d) in presence of the slab and aperture. A well - known 
phenomenon explained on the basis of the driven oscillator model [2]. Also, both vortices 
and a saddle point are exhibited by < S(r) > in the standing wave pattern between 
the nanocylinder and the slab. The vertical separation between the cylinder and the slit 
exit has been chosen in order to optimize this field enhancement. Figure 2(b) shows this 
field concentration in terms of the energy fiow magnitude | < S > | as the nanocylinder 
is gradually moved away from the slit at A = 366nm. In order to quantify this, we 
have averaged this quantity in an annulus surrounding the cylinder, whose internal circle 
coincides with the cylinder section. The image of Fig. 2(a) corresponds to the third peak 
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Figure 2: (a) Magnetic field modulus |H(r)| (colors in A/m units) and time aver- 
aged energy flow < S(r) >, (arrows in J/{m? ■ s units), maximum arrow length 
^ 80.02i^ey/(nm^ ■ s), minimum arrow length {)eV/{nm? ■ s)), localized on the surface 
of an Ag cylinder (radius R = 200nm, refractive index n = 0.186 + zl.61) in front of a slit 
in a slab at the same illumination as in Fig. 1(a). The distance between the cylinder 
surface and the exit plane of the slit is die = 3X/2 = 549?T,m (A = 366?T,m); (b) Variation 
of the concentration of | < S(r) > |, (in eV/lnm^ ■ s units), on the cylinder surface vs 
its vertical separation (in nm) from the slit exit, illuminated at A = 366nm (near the 
plasmon LSP51 resonance); (c) The same quantity vs wavelength at the distance seen 
in Fig. 2(a), die = 549nm. The black and green curves stand for the response of the 
slab alone and the red and blue ones for the slab with cylinder, respectively, (d) Electric 
field E(r) (in F/m units); both its norm (colors) and directions (arrows), under the same 
conditions as in Fig. 2(a). Black and red curves in Figs. 2(b) and 2(c) have been calcu- 
lated by averaging the quantity in an annulus of area A = 7r(((9/8)r)^ — r^) (see the two 
concentric circles in (a)) whose internal circle coincides with either the cylinder section, 
(case of the slit with the particle in Figs. 2(b) and 2(c)) or an imaginary circle coincident 
with that cylinder section, (case of the slit alone in Fig. 2(c)). Green and blue curves 
in Fig. 2(c) have been calculated by averaging the quantity in a rectangular monitor of 
area S = TGnmxbSAGnm = 4442. 96nm^ at the slit exit, (see the rectangle drawn in Fig. 
2(a)). 
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of Fig. 2(b) from the left. The spectrum of this enhancement vs the wavelength with the 
cylinder aX d = 549nm, can be seen in Fig. 2(c), which shows a comparison between this 
enhancement with the cylinder at the aforementioned distance d = 549nm with the one 
obtained in the same region with the slit alone (red and black curves, respectively), as well 
as in the region immediately outside the slit exit (blue and green curves). In addition. Fig. 
2(d) shows the pattern of E(r) for the same configuration as Fig. 2(a). This illustrates 
the enhancement of the electric field on the particle surface, as well as the interesting 
vortices described by its wavevector around it due to multiple reflections with the upper 
surface of the slab. This is in contrast with the case in which the particle is dielectric 
and the excited MDR is a whispering gallery mode (WGM), in which case the field is 
mainly confined inside the particle, and exponentially decays outside [ID]. In this latter 
case, the cylinder reflectivity is much lower and the multiple reflections are much weaker. 
This high reflectivity of the plasmonic cylinder produces not only a high concentration of 
transmitted energy around its surface, but also outside it. In particular, at the slit exit, 
(compare Fig. 1(a) and Fig. 2(a)). This is related to a relatively lower reflectivity at the 
slit entrance, and in fact, already this lower reflectivity constitutes another signature of 
a larger transmission into the space where the particle is. Hence, this cylinder acts as an 
extractor of transmittance upon the slit (see Figs. 2(a) and 2(d)). 




Figure 3: (a) Variation of the magnetic field concentration |H(r)| {in A/m units) on the 
cylinder surface vs. its horizontal separation (in nm) from the slit exit. The particle is 
illuminated at A = 366nm (i.e. near the nanocylinder LSP^i resonance). The black, red, 
green and blue curves stand for vertical distances between the particle surface and the slab 
equal to A/8, A/4, 3A/8 and A/2, respectively, (b) Map of |H(r)| (color) and < S(r) >) 
{J/{im? ■ s)) (arrows) when the cylinder is at a vertical and horizontal distance from the 
slit exit: A/8 and 17A/8, respectively, and illuminated at A = 364. 7nm (p - polarization, 
near the LSP51). Maximum arrow length ^ 26. 09 KeV/{nm'^ ■ s), minimum arrow length 
~ OeV/{n'm'^ ■ s). Calculations in Fig. 3(a) have been made like in Figs. 2(b) and 2(c). 

Figure 3(a) shows the behavior of both the magnetic field and the energy flow with a 
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lateral separation of the particle from the slit exit. Now, this asymmetry of its location 
produces that the standing LSP pattern on the cylinder surface disappears in its upper 
side, showing the energy circulation around its surface. Also, a standing wave is formed 
between the aperture and the reflecting cylinder which projects energy into the upper 
space. Again, the transmittivity of the slit is now much larger than when it is alone. 

2.3 Other sets of plasmonic nanocylinders 
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Figure 4: (a) Magnetic field modulus |H(r)| (in A/m) spatial distribution for five Ag 
cylinders, (radius R = 200?7,m, refractive index n = 0.186 + il.610), disposed in a bifur- 
cated chain, (bifurcation angle 6 = 45°, distance between cylinder surfaces dec = lOOnm), 
in front of the slit of Figs. 1(a) - (d) illuminated at A = SGA.Ynm (p - polarization, near 
the LSP51 resonance). The vertical distance between the bottom cylinder surface and 
the exit plane of the slit is Aq, (Aq = 366nm near the LSP^^i resonance of an isolated 
cylinder); (b) Time averaged energy flow < S(r) > (arrows in J/{m? ■ s units), maximum 
arrow length ^ Ib.QQKeV / {nm^ ■ s), minimum arrow length ^ QeV/ {nm? ■ s)) in the same 
configuration as in Fig. 4(a); (c) Electric field E(r) distribution, both norm (colors in 
V/m units) and directions (arrows) are displayed in the same configuration as in Fig. 
4(a). 

In order to test light transmission and concentration by means of aperture MDR - particle 



8 



LSP interaction and through adjacent cyhnder LSP couphng, Figs. 4(a) - (c) show |H(r)|, 
< S(r) > and E(r) for a bifurcated chain of Ag cyhnders near a sht. The distance between 
the sht and the cyhnders has been chosen so as to optimize the extraction of energy by 
the particles through the aperture. Light transport through the chain up to the upper 
particles is now obtained as a concentration of energy spread in the region of the four 
upper cylinders as shown in Figs. 4(a) - (c) with, of course, a lower intensity concentration 
in the aperture exit zone, (compare | < S(r) > | = 766.86ey/ (ram^ ■ s) in the area inside a 
rectangular monitor appropriately scaled and equivalent to that shown in Fig. 2(a)) over 
the slit exit when no particles are present against | < S(r) > | = 10698. 50eV^/(nm^ ■ s) 
averaged over the same monitor in the configuration of Figs. 4(a) - (c)). Notice that 
now, as before, the optimum resonant wavelength A = 364. 7nm for this combined system 
of particles and slit is again red - shifted with respect to the individual elements (A = 
359.4?2m). The field transmission through the particle chains is accomplished by both 
the propagating waves surrounding the set and a field hopping process between neighbor 
particles. An appropriate choice of set parameters and illumination, allows one to select 
the particle of the set with the most enhanced LSP field intensity on its surface. 

Chains of smaller size nano cylinders at distance dec from each other, show dipolar plasmon 
coupling between neighbors, [see |H(r)| in Figs. 5(a) and 5(b)], both near the exit plane 
of the slit and at a larger distance from it. A detail of the transmission in the first 
three cylinders of such a chain can be seen in Figs. 5(c) - (d), where the energy flow 
and the electric field are shown. The latter quantity exhibits on the surface of each 
cylinder a typical dipolar LSP distribution, [see Fig. 5(d)]. This transmitted light through 
the particles does not take place by tunneling, but by a hopping mechanism, like in 
dielectric molecular photonic crystal rows [TTj, due to a dipolar interaction (where dec 
is the distance between cylinder surfaces) with frequency splitting of the single plasmonic 
cylinder spectral line, and a redshift of its extinction peak. This coupling diminishes the 
energy of the ensemble in the bonding state for this configuration (parallel dipoles) [5]. 
The field confinement along the chain line and between cylinders in the case in which 
they are close to the slit [cf. Figs. 5(a), (c) and (d)] contrasts with the large diffraction 
occurring when their distance to the slab is large [cf. Fig. 5(b)], this latter case also 
showing a strong standing wave below the chain. The energy flow | < S > | in the area 
immediately outside the slit exit is reduced with respect to the case of the slit alone in the 
case of Fig. 5(a) (compare for Fig. 5(a) 192.48eV^/(nm^-s) against 1536.35ey/(nm^-s) for 
the same slit alone) however, it is enhanced in the case of Fig. 5(b) (24674. 90el^/ (nm^ ■ s) 
compared with 6531. 39eV/ (rim^ ■ s) for the same slit alone. These numbers were obtained 
by averaging in a suitable rectangular monitor at the slit exit. Nevertheless, both linear 
configurations render a good response regarding to the averaged energy flow concentrated 
around the cylinders. (See the corresponding low reflected energy in Fig. 5(b)). The 
linear chain near the slit of Fig. 5(a) achieves an averaged | < S > | of llOOeV / {nrri^ ■ s) 
against 188.876^/ (nm^ ■ s) in the same area when the slit transmits without the cylinders. 
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Figure 5: (a) Magnetic field modulus |H(r)| [A/m) distribution for five Ag cylinders 
(radius R = 30nm, refractive index n = 0.173 + 21.95) of a linear chain, (distance between 
cylinder surfaces dec = lOOnm), in front of a slit in a W slab, (slit width d = 39.59?T,m, 
slab width D = 2610?7,m, slab thickness h = 237.55nm, refractive index n = 3.39 + z2.41) 
illuminated at A = 400nm (p - polarization). The vertical distance between the bottom 
cylinder surface and the exit plane of the slit is Ao/8 ( Aq = 346nm is a wavelength near 
the LSPii of an isolated cylinder); (b) The same configuration as in Fig. 4(a), changing 
the distance from the slit to 3Ao/2 (Aq = 346nm), illuminated at A = 349. 3nm (p - 
polarization); (c) Detail of the time averaged energy flow < S(r) > (in J/m^ ■ s) in the 
first three cylinders of the chain for the arrangement of Fig. 4(a), both the norm (colors) 
and directions are shown, (maximum arrow length ^ 1.7 5 KeV/{nm'^ ■ s), minimum arrow 
length ^ OeV/{nm'^ ■ s)). (d) Detail of the electric field E(r) in the first three cylinders of 
the chain for the case of Fig. 4(a), both the norm (colors in V/m) and directions (arrows) 
are shown. 
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Figure 6: Magnetic field modulus |H(r)| (colors in A/m units) and time averaged energy 
flow < S(r) > (arrows in J/m?-s units, maximum arrow length ^ 1.52-10^i^ey/(nm^-s), 
minimum arrow length OeV/ {nim?-s)) concentrated on the surfaces of seven Ag cylinders 
(radius R = 30nm, refractive index n = 0.186 + zl.61) in bifurcation, (bifurcation angle 
9i = 45°, distance between cylinder surfaces dec = Onm), with elbow, (elbow angle 02 = 
45°), in front of a slit in a W slab (slit width d = 39.59?7,m, slab width D = 2Q10nm, slab 
thickness h = 237. 55nm, refractive index n = 3.39 + z2.66), illuminated at A = 364. Tram 
(p-polarization) which is near the LSPu of each cylinder. The vertical distance between 
the bottom cylinder surface and the exit plane of the slit is die = Aq = 346nm, near the 
LSPu of the isolated cylinder. 
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Other distributions considered include those of bifurcations with elbows of nanoparticles 
at subwavelength distance dec from each other, as illustrated in Figure 6, which shows 
(|H(r)| in color and < S(r) > in arrows. Now the dipolar interaction is of order 
since it takes place in the near field. The enhancement of transmitted light on top of 
the set is now quite sharp as the cylinders approach each other (not shown). This set 
presents the bottom cylinder placed at a distance from the slit that optimizes the energy 
concentration both on the top particles and at the aperture exit. Most light passes through 
the aperture and reaches the cylinder set, although the radiation pattern reveals coupling 
between cylinders and slab. Notice the low value of |H(r)| in the reflection region below 
the slab in Fig. 6. The standing wave between the slab and the cylinders now embraces 
the set in one of the maxima, and makes it to strongly emit upwards along three main 
directions, hence exhibiting a nanoantenna-like behavior. (Radiation directions can be 
controlled by changing the configuration at the chosen wavelength). Again, this effect 
is more pronounced with the cylinders at a certain distance from the slit, as shown in 
Fig. 6, than when they are very close to it. In either case, however, this intensity, 
strongly transmits even into the aperture exit [42216. 07 eV/{nm'^ ■ s) for Fig. 6 against 
2783.01e1//(nm2 ■ s) for the same slit alone), (these values obtained by an average over 
the area of a suitable rectangular monitor at the slit exit, although most energy is now 
transmitted to the cylinders, and associated with a very small reflection below the slab. 
This process of transmission up to the top of the set is much more efficient than that of 
transmission by the excitation of the morphological resonance of the slit alone, which was 
shown in Fig. 2(c). 

2.4 A metallic photonic crystal in front of an array of slits 

Aperture supertransmission is further enhanced by periodically repeating it in the slab [10] 
- [13]. Accordingly, we next arrange chains of metallic particles, each placed in front of a 
slit of such an array. Introducing some distance between cylinders in each chain, we build 
in this form a metallic photonic crystal (PC). 

Figures 7(a) and 7(b), corresponding to the PC alone and to the PC in front of the 
slit array, respectively, show the first step of our approach. At the chosen illumination 
wavelength A = 400nm the slit grating is supertransmitting and so does the array because 
this wavelength is both close to a Rayleigh resonance of the slit array and near the LSPu 
of the PC cylinders. As a result, a strong concentration of field in the vertical rows of 
the PC appears. This is again due to a dipolar interaction between neighbor particles of 
the same vertical row and between adjacent rows as discussed in connection with Figs. 
5(a) - (d). The enhancement of the field is alternating between particle gaps in each row, 
at difference with the case of a single chain, shown in Figs. 5(a) - (d). The qualitative 
aspect of the distribution of Figs. 7(a) and 7(b) is similar, except for the coUimation effect 
produced by the slit grating, and the slight weakening of the energy enhancement between 
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Figure 7: (a) Magnetic field modulus |H(r)| [A/m) distribution in a photonic crystal 
(horizontal period = 275.50nm, vertical period ay = 75nm) formed with 68 Ag cylin- 
ders (radius R = 30nm, refractive index n = 0.173 + zl.95) and illuminated at A = 400nm 
(p - polarization, the LSPn excited is near that of the isolated particle which would ap- 
pear at A ~ 346nm); (b) Magnetic field modulus |H(r)| distribution in the Ag photonic 
crystal of Fig. 7(a), this time in front of a grating of slits practiced in a slab (period 
P = 275.50nm, slit width d = 39.36nm, slab width D = 8P, slab thickness h = dAAQnm, 
refractive index n = 3.39 + i2.41), illuminated at A = 400nm (i.e. near the LSPu of an 
isolated particle), (p - polarization). The distance between the cylinder surfaces of the 
first row and the exit plane of the slits is 22.5nm (fitted to get the best response). 
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particles due to intensity concentration within the shts, shown in Fig. 7(b). These results 
in Fig. 7(b) contrast with those in which the incident wavelength is out of resonance of 
the cylinders, then light passes through the PC with scarce interaction with the cylinders 
(not shown here). 
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Figure 8: (a) Norm | < S > | in eV/{nm? ■ s) vs. wavelength A (in nm) transmitted 
by the W array of slits alone of Fig. 7(b). (b) Average energy flow norm | < S > | 
in eV/ {nm?' ■ s) units vs. wavelength passing through the PC of Ag cylinders of radius 
R = 30nm, [horizontal period = 275.50?7,m, vertical period ay = IQOnm, see also Figs. 
9(a) - (c)]. (c) Norm (| < S > |) vs. wavelength for the same PC, now in front of the 
W grating of Figure 8(a). The black squared and red circled curves stand for the norm 
of energy flow averaged over each square (120 x 120?T.m^) circumscribed to each cylind.er 
section and over each rectangular strip (120 x 920nm^) circumscribed to each PC vertical 
row, respectively. In the case of Figure 8(a), these circles are imaginary and coincide with 
the cylinders sections of Figures 8(b) and 8(c), (see also Figs. 9(a) - (c)). 

In the second step of our analysis, the distance between cylinders in each vertical chain 
is incremented to a value comparable to that of the horizontal distance between chains. 
Figures 8(a) - (c) correspond to the responses of the grating alone, this new PC alone 
and the combination of both, respectively. Figure 8(b) shows an effective bandgap in the 
XM direction of transmission in PC reciprocal space (upwards direction in Figs. 7(a) 



14 



- (b) and Figs. 9(a) - (c)). Nevertheless, an enhancement of transmission along this 
direction is achieved in the whole range studied (compare the values of Fig. 8(b) to those 
of Fig. 8(c)). Furthermore, a transmission peak rises for the PC/grating arrangement, 
in the range where the gap was for the PC alone (see Figs. 8(c) and 8(b)), the two last 
effects being due to the presence of the grating. Note the match in wavelength between 
the transmission peak of the grating alone in Fig. 8(a) and that mentioned above for the 
combination of PC and grating in Fig. 8(c). 




Figure 9: (a) Map of the magnetic field norm |H(r)| (in A/m units) for the Ag cylinder 
photonic crystal, studied in Fig. 8(b), placed in front of the exit of the W array of Figs. 
7(b) and 8(a), illuminated at A = AbOnm (see Fig. 8(c)). The distance between the 
first horizontal row lower edges and the slab is 22.5nm. (b) Detail of the time averaged 
energy flow S(r) (in J /im?-s units), showing both its norm (color) and directions (arrows), 
maximum arrow length 11. 23 KeV/{n'm'^ ■ s), minimum arrow length ^ QeV/{nm? ■ s)). 
(c) Detail of the electric field E(r) (in V/m units), both norm (color) and directions 
(arrows) are shown. 

By contrast to Fig. 7(b), when the slit grating is added to this PC, and illumination 
is out of both the LSP resonance and the PC effective gap, (see Fig. 8(b)), one again 
observes an enhancement of transmission into the PC, a strong intensity concentration in 
the slits, and a relatively low reflection below the slab (as example, that shown in Figs. 
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9(a) - (c)), even though now this reflection is larger, relative to the transmittivity of the 
slits, than that shown for resonantly excited nanoparticle sets in front of one slit (compare 
with Figs. 5(a) - (d) and 6). On the other hand, the presence of the slits modifles the gap 
of the PC which, as shown in Fig. 8(c) has a transmission peak near A = 360?7,m. This is 
a new feature of the effect of the slits upon the transmittivity of the PC. 



16 



3 Conclusion 

We have presented a theoretical and simulation study which shows how the excitation 
of localized surface plasmons of nanoparticles in front of subwavelength apertures, {Ag 
cylinders and W slabs have been used) under polarization that excites the aperture prop- 
agating modes (p - polarization in 2D), may enhance the transmission. Sets of plasmonic 
cylinders exhibit coupling and propagation through their elements. There is a specific 
distance for each of such cylinder - slit configurations that optimizes the transmitted en- 
ergy passing into the particles. This suggests a control of light transmission via plasmonic 
particles. Also, we observe that it is possible to fitting particle set parameters and illu- 
mination in such a way that the transmitted intensity is concentrated in certain cylinders 
when the stationary regime of propagation has been reached. 

The case of a metallic photonic crystal {Ag PC) in front of a metallic slit array shows the 
effects of effective bandgaps on the interparticle LSP transmission, and the enhancement 
of the transmittivity of the grating due to the excitation of LSPs in the metallic PC. All 
these results are also expected with 3D particles in front of apertures with any geometry, 
in particular subwavelength holes, and have a potential for controlling transmitted near 
fields at the nanoscale. 
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